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INTRODUCTION 


Hydrometeor  concentrations  aloft  are  determined  by  radar 
and  by  airplane-mounted  samplers.  Radar  measurements  are  stronqly 
biased  toward  larger  particle  sizes,  and  do  not  discriminate  between 
hydrometeor  types.  Airplane  samples  can  be  used  to  calibrate  the 
radar  measurements,  identify  hydrometeor  types,  and  fill  out  the 
small-particle  tails  of  the  size  distributions. 

One  of  the  problems  with  airplane  sampling  of  hydrometeors 
is  that  hydrometeor  concentrations  tend  to  be  distorted  at  sampling 
points  near  the  fuselage  surras*.  This  report  is  concerned  with 
theoretical  assessment  of  this  distortion  for  specific  airplane- 
mounted  instruments.  Theoretical  bases  of  the  method  used  and 
several  applications  are  presented  in  refs.  1 and  2.  A summary  of 
the  theory  is  given  in  the  next  section. 

In  this  report,  additional  support  for  use  of  steady- 
state  drag  coefficients  for  calculation  of  accelerative  particle 
motion  is  discussed.  Assessments  of  sampling  adequacies  of  several 
instuments  mounted  on  a Learjet  36  airplane  are  reported.  This 
includes  assessment  of  a total  water  content  indicator  in  the  nose 
of  the  Learjet  that  requires  computation  of  flow  into  the  airplane 
nose:  a new  application  of  the  method.  Finally,  results  of  a pre- 
liminary attempt  to  prepare  a generalized  graph  to  be  used  to  esti- 
mate "safe"  sampling  distances  for  fuselage-mounted  instruments  are 
compared  with  results  obtained  from  a number  of  specific  studies. 


1.  H.Gi  Norment  and  R.G.  Zalosh,  "Effects  of  Airplane  Flowfields  on 
Hydrometeor  Concentration  Measurements",  AFCRL-TR-74-0602  (6  Dec. 
1974).  AD-A006  690. 

2.  H.G.  Norment,  "Effects  of  Airplane  Flowfields  on  Cloud  Water  Con- 
tent Measurements",  AFCRL-TR-75-0231  (30  April  1975).  AD-A014  807. 
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BACKGROUND 
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NATURE  OF  THE  PROBLEM 

Particle  sampling  devices  of  interest  here  usually  measure 
particle  flux  through  a small  space  adjacent  to  an  airplane  fuselaqe. 
Unless  the  sampling  space  is  distant  enough  from  the  airplane  to  be 
in  the  free-stream,  particles  of  certain  sizes  will  interact  with  the 
flow  about  the  fuselage  to  cause  flux  distortion.  The  situation  is 
illustrated  in  Fig.  1 for  100  urn  diameter  water  drops  in  airflow  about 
a prolate  ellipsoid  of  fineness  ratio  2.  Note  the  impaction  on  the 
ellipsoid  of  the  drop  closest  to  the  ellipsoid  symmetry  axis,  and  note 
the  substantial  deflections  of  the  next  closest  trajectories.  Tra- 
jectory deflection  causes  high  particle  concentrations  and  high  con- 
centration gradients  to  be  observed  at  a point  such  as  the  one  marked 
01  in  the  figure.  At  point02,  deflection  and  impaction  combine 
to  produce  a region  void  of  particles,  a so-called  "shadow  zone". 
Smaller  drops,  with  much  less  inertia,  tend  to  follow  the  airflow 
more  exactly  such  that  less  distortion  is  observed.  Drops  larqe 
enough  to  have  very  high  inertia  substantially  ignore  the  airflow, 
and  again  little  distortion  is  observed.  Therefore,  distortion  is 
significant  over  a limited  range  of  particle  sizes. 

CONCENTRATION  FACTOR 

Principal  results  of  this  work  are  expressed  in  a quantity 
called  concentration  factor.  Concentration  factor,  Cp,  is  defined  as 
the  ratio  of  particle  flux  at  the  sampling  or  target  point,  F^,  to 
the  particle  flux  in  the  freestream,  F, 
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FLOW  DIRECTION 


Trajectories  of  100  pm  diameter  water  drops  in  potential  airflow  about  a prolate  ellipsoid 
fineness  ratio  2.  (The  ordinate  scale  is  expanded  by  a factor  of  2.) 


A 


The  ratio  of  the  particle  concentration  at  the  target  point  to  free- 
stream  concentration,  C^,  is 

CM  > CF  V/Vt  • (2) 

where  V is  free-stream  airspeed  and  is  airspeed  at  the  tarqet 
poi nt.  Usual  ly  V/V^.  - 1. 

Concentration  factor  is  determined  by  computing  particle 
trajectories  from  the  freestream  (initial  plane  in  Fig.  2)  to  a 
small  area  in  the  target  plane  (Fig.  2)  that  surrounds  the  target 
point  such  as  to  define  a particle  flux  tube.  Since  the  particle 
mass  transfer  rate  through  the  tube  is  constant  at  all  cross  sec- 
tions, it  is  easily  shown  that 


where  A and  are  the  flux  tube  cross  sectional  areas  in  the  free- 
stream and  at  the  target  point.  In  the  limit  as  A and  A^  approach 
zero,  eq.  (3)  becomes  exact. 

PARTICLE  TRAJECTORY  CALCULATION 

We  assume  that  the  bulk  air  flow  is  not  perturbed  by  the 
particles.  Moreover,  since  particle  density  is  large  compared  to 
that  of  air,  we  can  neglect  buoyancy  and  inertial  reaction  of  the 
fluid  to  obtain  the  three-dimensional,  normalized  equation 


^ i j i_  (?  ) !mVs  _ t 

fh  Ivs  1 8 BN,sRH 


w 


a 


T 


still-air  terminal  settlina  speed 

of  the  particle 

unit  vector  in  the  z (upward) 

direction 

time 


Fn  * V2/ (Lg) 


Froude  Number 


Rw  = — v -v  I V 
N n a p 


BN  = cdV 


Dimensioned  quantities  are: 


Reynolds  Number 

Davies  Number 

particle  drag  coefficient 


6 

P 

n 

g 

v 

L 


particle  dimension 
air  density 
air  viscosity 

gravity  acceleration  constant 
freestream  airspeed 
a characteristic  dimension  of  the 
fuselage 


Here  length  is  normalized  by  L,  velocity  by  V and  time  by  L/V.  RN 

• 1 > s 

and  are  for  still-air,  terminal  settling  of  the  particles. 

Starting  at  the  initial  plane,  eq.  (4)  is  integrated 
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with  respect  to  time  in  three-dimensional  space  via  the  code  DVDP 
(3) 

of  Kroglv  ' till  the  target  plane  is  reached.  The  method  used  to 
compute  v,  at  each  time  step  is  described  next,  and  then  drag  coef- 

a 

ficients  are  discussed. 

THREE-DIMENSIONAL  FLOW  CALCULATION 

In  performing  concentration  factor  calculations  for  sam- 
pling sites  on  particular  airplanes,  it  is  important  to  use  three- 
dimensional  airflow.  This  is  the  only  way  to  properly  account  for: 
airplane  oeometry  and  angle-of-attack,  airspeed,  altitude,  and  par- 
ticle settling. 

Cloud  physics  airplanes  are  subsonic,  sampling  runs  being 
made  typically  between  100-150  kts  indicated  airspeed  (IAS).  Par- 
ticle measurement  points  are  beyond  the  skin- friction  boundary  layer, 
and  should  be  placed  to  avoid  regions  of  separated  flow.  Therefore, 
potential  (i.e.  , frictionless,  incompressible,  laminar)  flow  calcula- 
tions are  quite  adequate.  We  use  a code  developed  by  Hess  and  Smith/4, 5 
for  calculating  potential  flow  about  arbitrary  three-dimensional 
bodies.  The  Hess-Smith  code  requires  input  of  a digital  description 
of  the  fuselage  surface.  This  consists  of  the  coordinates  of  the 
corner  points  of  a large  number  of  contiguous,  plane,  quadrilaterals. 

An  example  of  the  digital  description  of  a fuselage  is  shown  in  Fig.  3. 

Accuracy  of  the  flow  calculations  has  been  checked  with 
excellent  results  by  Hess  and  Smith/4)  for  many  bodies  for  which  ana- 
lytical solutions  are  available.  Norment  and  Zalosh/*)  have  compared 
computed  trajectories  around  ellipsoids  in  analytical  flow  fields  with 
similar  trajectories  in  Hess-Smith  flow  fields,  and  they  have  compared 
current  trajectory  results  with  prior  work;  agreement  is  excellent. 

3.  F.T.  Krogh,  "Variable  Order  Integrators  for  Numerical  Solutions  of 
Ordinary  Differential  Equations,"  Jet  Propulsion  Lab  Technoloay 
Utilization  Document  No.  CP-2308  (November  1970). 

4.  J.L.  Hess  and  A.M.O.  Smith,  "Calculation  of  Non-Lifting  Potential 

Flow  About  Arbitrary  Three-Dimensional  Bodies,"  McDonnel  Douglas 
Report  E.S.  40622  (15  March  1962).  AD-282  255. 

5.  J.L.  Hess  and  A.M.O.  Smith,  "Calculation  of  Potential  Flow  About 
Arbitrary  Bodies,"  Progress  in  Aeronautical  Sciences,  Vol . 8,  edited 
by  D.  Kuchemann  (Pergammon  Press,  New  York,  1967). 
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(b) 


Figure  3.  Computer  prepared  plots  of  the  digital  description 
of  the  Learjet  36  fuselage:  (a)  Forward  fuselage,  (b)  ComDlete 
fuselage.  Instrument  locations  are:  ASP#  , Replicator#  , 
Snow  Stick  ▼ . 
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PARTICLE  DRAG  COEFFICIENT 


Davies^shows  that  still-air  terminal  settlinq  of  spheres 
can  be  generalized  -in  terms  of  the  dimensionless  numbers  R()  c and 
..  Over  the  range  from  the  smallest  spheres,  which  settle  under 
viscous  flow  conditions  and  obey  Stokes  law,  to  spheres  much  larger 
than  of  interest  here,  and  for  any  Newtonian  fluid,  a reproducible 
single-valued  relationship  between  R,.  and  B,,  exists.  Furthermore, 

<1  , S li  y S 

B.j  s is  independent  of  settling  speed,  being  a function  of  fluid  and 
sphere  properties  only;  thus  for  given  sphere  and  fluid,  Pfj  _ and 
hence  V can  be  calculated.  Polynomials  by  which  R..  can  be  com- 

^ 1 y S 

puted  as  a function  of  B„  were  derived  by  Davies  from  a composite 

IN , S 

of  many  sets  of  experimental  data. 

Since  the  work  of  Davies,  it  has  been  found  repeatedly  that 
this  treatment  is  applicable  to  particles  of  other  shapes,  providing 
settling  is  steady  and  particle  orientation  is  stable. 

For  the  trajectory  calculations  required  here,  the  problem 
must  be  turned  around.  In  addition  to  gravity  settling,  there  is  a 
particle  velocity  component  (relative  to  air)  caused  by  the  distur- 
bance of  the  passing  airplane.  At  any  time  step  in  the  numerical 
integration  of  eq.  (4),  v^-v^  (and  hence  R^)  is  known,  and  must  be 
determined.  For  viscous  motion  (i.e.  , Stokes  flow,  where  R^  < 1) 

B^  = 24  R^  and  eq.  (4)  can  be  integrated  without  question.  However, 
for  larger  R^  the  steady-state  drag  data  determined  experimentally  for 
terminal  settling  must  be  used  to  compute  accelerative  particle  motion. 

Experimental  measurements  by  Keirr/^  and  a theoretical  analy- 
sis by  Crowe,  et  al . ^ indicate  that  if  the  acceleration  modulus. 


6.  C.N.  Davies,  "Definitive  Equations  for  the  Fluid  Resistance  of 
Spheres,"  Proc.  Phys.  Soc.  (London)  57,  259-270  (1945). 

7.  S.R.  Keim,  "Fluid  Resistance  to  Cylinders  in  Accelerated  Motion," 

J.  Hydraulics  Div. , Proc.  Amer.  Soc.  Civil  Enq. , 6,  oaoer  1113  (1956). 

8.  C.T.  Crowe,  J.A.Nicholls  and  R.B.  Morrison,  "Drag  Coefficients  of 
Inert  and  Burning  Particles  Accelerating  in  Gas  Streams,"  Ninth  Symp. 
(Int'l.)  on  Combustion,  Academic  Press,  pp.  395-405  (1963). 
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steady-state  drag  coefficients  can  be  used 

without  significant  error  to  compute  accelerative  motion.  has 

never  been  found  to  exceed  10- : in  our  trajectory  calculations. 

For  small  water  drops,  which  are  spherical,  the  polynomial 

eguations  of  Davies^  are  used  to  compute  V.  while  for  larger,  dis- 

s (91 

torted  drops,  the  equations  of  Foote  and  du  Toit'  ' are  used.  To  com- 
pute B.,  from  Rn,  inverse  polynomials  have  been  developed,  using  the 

N (101 
data  set  given  by  Davies  for  small  drops  and  Gunn  and  Kinser'  for 

large  drops. 

Concentration  factors  also  are  calculated  for  a variety  of 
ice  crystal  forms  and  for  broad  ranges  of  ice  crystal  sizes.  Ice 
crystal  drag  and  dimensional  data  are  described  in  refs.  1 and  2. 


AN  = 4 


dV 

dt 


/VP 


is  smaller  than  about  10~2, 


9.  G.B.  Foote  and  P.S.  du  Toit,  "Terminal  Velocity  of  Raindrops  Aloft," 
J.  Appl . Meteor.  8,  249-253  (1969). 

10.  R.  Gunn  and  G.D.  Kinser,  "The  Terminal  Velocity  of  Fall  for  Water 
Droplets  in  Stagnant  Air,"  J.  Meteor.  6,  243-248  (1949). 
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LEARJET  36 


INSTRUMENT  LOCATIONS 

A Learjet  36  airplane  (Fig.  4)  has  been  instrumented  for 
cloud  physics  research  by  Meteorology  Research,  Inc.  (MRI)  as  shown 
in  Fig.  5.  Instruments  of  interest  here  are: 

Total  Water  Content  Indicator  (TWCI) 

Particle  Measuring  Systems  (PMS)  1-D  Particle  Spectrometers 

Snow  Stick 

Formvar  Replicator 

PMS  Axial  Scattering  Probe  (ASP) 

Engineering  drawings  of  surface  contours  of  the  fuselaae 
were  obtained  from  the  Gates-Lear jet  Corp.  The  engineering  coordi- 
nate system  is  as  follows  for  this  airplane: 

FS  Along  the  fuselage  axis,  increasing  in  the  aftward 
direction,  with  origin  86.75"  forward  of  the  nose  tip. 

BL  Horizontal  direction  perpendicular  to  the  fuselage 
axis  with  origin  at  the  axis. 

WL  Vertical  direction,  increasing  upward,  with  origin 
at  the  lowermost  level  of  the  fuselage. 

For  the  calculations  reported  here,  the  origin  is  shifted 
to  the  nose  tip  (actually  to  the  center  of  the  TWCI  orifice)  and  all 
coordinates  are  normalized  by  dividing  by  the  radius  of  the  upper 
portion  of  the  cylindrical  fuselage  section  (32").  Coordinates  of 
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the  instrument  sampling  points  are  given  in  Table  1. 

An  accurate  digital  description  of  the  fuselage  was  devel- 
oped from  the  engineering  drawings.  Fig.  3 shows  computer  plots  of 
the  fuselage,  which  is  represented  by  1104  quadri lateral s. 

Particle  trajectory  calculations  were  done  for  30  kft 
altitude,  for  which  flight  conditions  are: 


Indicated  airspeed 
True  airspeed 
Angle  of  attack 
Temperature 
Air  density 


165  kt  (84.9  m sec" 1 ) 
270  kt  (138.8  m sec-1 ) 
4.5  deg 
228.8°  K 
0.45904  kg  m"3 


RESULTS  FOR  THE  SNOW  STICK,  REPLICATOR  AND  ASP 

As  shown  in  Fig.  5,  these  instruments  are  mounted  through 
the  cabin  windows.  Distances  from  the  fuselage  skin  are  13  or  14 
inches,  and  considering  the  slenderness  of  the  fuselage  (radius  32  in.) 
are  well  clear  of  shadow  zones.  Table  2 shows  concentration  factors 
as  a function  of  water  drop  diameter  at  the  instrument  sampling  loca- 
tions. To  bracket  the  shadow  zones,  concentration  factors  were  com- 
puted at  points  successively  closer  (by  2.54  inch  increments)  to  the 
fuselage  relative  to  the  instrument  sampling  points.  Results  are 
shown  in  Table  3. 

Consistent  with  being  the  furthest  aft,  the  snow  stick  has 
the  highest  concentration  factors  and  the  deepest  shadow  zone.  To 
determine  whether  or  not  the  snow  stick  is  shaded  for  ice,  concentra- 
tion factors  were  computed  for  solid  and  hollow  ice  columns  and  for 
hexagonal  plates.  Tables  4 and  5 give  the  ice  crystal  properties 
and  Table  6 gives  the  concentration  factor  results.  No  shadowing  is 
indicated. 
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TABLE  2 


WATER  DROP  CONCENTRATION  FACTORS 
AT  TiiE  LEARJET  36  SNOW  STICK,  FORMVAR  REPLICATOR  and  ASP 


Water 

Drop 

Diameter 

(iLm) 

50 
80 
100 
300 
son 


CONCENTRATION  FACTOR,  Cp 


Snow  Stick 

Repl icator 

ASP 

1.148 

1.155 

1.151 

1.203 

1. 173 

1.136 

1.  194 

1. 153 

1.116 

1.046 

1.032 

1.023 

0.987 

0.988 

0.998 

17 


PROPERTIES  OF  COLUMNAR  ICE  CRYSTALS 


18 


TABLE  5 

PROPERTIES  OF  HEXAGONAL  PLATE  (Pla)  ICE  CRYSTALS 


Diameter  of 


Diameter 

(ym) 

Thickness 

(um) 

Mass 

(pq) 

Water  Drop  of 
Equal  Mass 
(’Pm) 

Density 

(kq/m3) 

10 

5.68 

% 2.582  x 10"4 

% 7.9 

* 700. 

30 

9.30 

o 3.591  x 10* 3 

% 19 

•y  660.6 

50 

11.70 

% 9.203  x 10-3 

* 26 

*y  484.0 

100 

15.97 

% 3.351  x 10-2 

* 40 

% 323.0 

300 

26.16 

0.4779 

97 

312.5 

500 

32.90 

2.145 

160 

401.5 

1000 

44.91 

20.22 

338 

693.1 
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TABLE  6 


ICE  CONCENTRATION  FACTORS  AT  THE  LEARJET  36  SNOW  STICK. 


a.  COLUMNS 


Concentration 

Factor,  Cp 

Length 

(urn) 

Sol  id 

Hoi  low 

20 

1.032 

30 

1.048 

1.039 

50 

1.089 

1.060 

100 

1.181 

1. 116 

300 

1.241 

1.211 

500 

1.221 

1.242 

1000 

1.164 

1.251 

3000 

1.086 

1.178 

b.  HEXAGONAL  PLATES 


Diameter 

(um) 

CF 

30 

1.050 

50 

1.053 

100 

1.062 

300 

1.122 

500 

1.190 

1000 

1.262 
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PARTICLE  TRAJECTORY  ANGLE  AT  THE  1-D  PARTICLE  SPECTROMETERS 


J 

As  shown  in  Fig.  5,  a pair  of  PMS  particle  spectrometers 
are  mounted  above  the  fuselage  nose.  The  sampling  volumes  are  18" 
aft  of  the  nose  tip  and  are  sliqhtly  more  than  8"  from  the  nose  surface. 

We  expected  no  significant  concentration  distortion  and  found 
none,  but  there  was  interest  in  determining  trajectory  angles  relative 
to  the  horizontal  at  the  sampling  points.  Figure  6 shows  a short  sea- 
ment  of  trajectory  for  300  urn  diameter  water  drops  to  one  of  the  spec- 
trometers. The  airplane  axis  is  tilted  at  4.5°.  The  trajectory 
is  almost  horizontal,  being  tilted  in  the  same  sense  as  the  airplane 
axis  at  an  angle  less  than  0.5  at  the  sampling  point. 

TOTAL  WATER  CONTENT  INDICATOR 

Instrument  Description 

The  intake  orifice  of  the  TWCI  consists  of  a 2 inch  diameter 
hole  in  the  nose  tip  of  the  fuselage.  Intake  air  velocity  is  about 
30  ft  sec'1,  but  will  vary  somewhat  with  flight  conditions.  For  the 
calculations  reported  here,  we  have  used  an  intake  of  5 of  the  free- 

stream  airspeed,  which  corresponds  to  23  ft  sec'1  at  30  kft  altitude. 

o * 

A 6 angle  of  attack  was  used  for  the  TWCI  calculations. 

As  illustrated  by  Fig.  3 the  fuselage  is  represented  by  a 
large  number  (1104)  of  plane,  contiguous  quadrilaterals.  Usually  the 
fuselage  symmetry  is  used  so  that  only  half  of  the  total  number  of 
quadrilaterals  need  be  specified.  However,  this  requires  that  quad- 
rilateral boundaries  be  defined  along  the  symmetry  plane,  and  unfor- 
tunately this  plane  bisects  the  TWCI  orifice.  Since  flow  accuracy  is 
highest  at  quadrilateral  centers  and  lowest  at  boundaries  we  consid- 
ered it  unwise  to  use  the  fuselage  symmetry  and  have  quadrilateral 
* 

At  the  time  the  TWCI  calculations  were  done,  the  best  information 
available  indicated  a 6 A0A,  whereas  a value  of  4.5°  should  have 
been  used  to  be  consistent  with  the  other  calculations  reported 
here.  However,  we  find  that  the  A0A  discrepancy  produces  negligi- 
ble differences  in  the  Cr  results. 
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Figure  6.  Computer  prepared  plot  showinq  a seqment  of 
trajectory  for  a 300  um  diameter  water  drop  to  one  of 
the  PMS  1-D  particle  spectrometers.  The  airplane  axis 
is  inclined  at  4.5  degrees  to  the  horizontal. 
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boundaries  running  through  the  orifice  center.  Thus  a special  ver- 
sion of  the  nose  was  developed  for  the  TWCI  study  which  does  not 
use  the  symmetry.  This  version,  shown  in  Fig  7a,  allows  the  center 
of  the  orifice  to  coincide  with  a quadrilateral  center.  Since  not 
using  the  symmetry  essentially  doubles  the  time  for  each  air  veloc- 
ity computation,  the  fuselage  is  described  faithfully  only  as  far 
aft  as  the  base  of  the  cockpit  windshield.  Beyond  this  point  the 
fuselage  cross  section  is  kept  constant  as  shown  in  Fig.  7a.  This 
description  of  the  nose  requires  471  quadrilaterals. 

Figure  7b  shows  an  enlargement  of  the  nose  tip  of  the 
Learjet  fuselage.  The  2 inch  diameter  orifice  is  represented  by 
a plane  consisting  of  21  quadrilaterals.  The  airflow  computer  code 
was  modified  such  that  each  of  these  21  quadrilaterals  "leaks" 
inward  a specified  portion  (in  this  case  5%)  of  the  freestream  flow. 

Flowfield  Studies 

Since  we  had  not  previously  computed  flow  into  and  around 
an  orifice,  we  performed  a detailed  study  of  the  flow  field  in  the 
vicinity  of  the  nose  tip.  Calculations  were  done  for  zero,  5°/  and 
100 ; intake  relative  to  freestream  airspeed.  Flowfields  in  the  fuse- 
lage symmetry  plane  (i.e. , plane  y=0)  are  shown  in  Fig.  8.  (Note 
that  there  is  little  apparent  difference  between  the  zero  intake  and 
5'  intake  cases.)  Careful  study  of  the  numerical  results  indicate 
that  the  intake  is  accurate  and  uniform  over  essentially  the  entire 
surface  of  the  orifice. 

Figure  9 shows  a stereo-view  of  100  urn  diameter  water  drops 
to  the  nose  with  5%  intake  into  the  orifice.  In  the  case  of  these 
particles  the  trajectories  deviate  little  from  straight  lines. 

Concentration  Factors,  Concentration  Ratios  and  Concentration 
Correction  Factors 

Concentration  factors  were  calculated  in  the  usual  manner 


k 
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Figure  7.  Special  digital  description  of  the  Learjet  36 
nose  for  the  TWCI  study:  (a)  Complete  nose  description 
(b)  Enlargement  of  the  nose  tip.  The  center  of  the  TWCI 
orifice  is  marked  • . 


Figure  8b.  Flowfield  around  the  Learjet  36  nose  tip  in  the  plane 
y=0  for  intake  of  5%  of  the  freestream  flow  into  the  TWCI  orifice 


////////i 


Figure  8c.  Flowfield  around  the  Learjet  36  nose  tip  in  the  plane 
y=0  for  intake  of  100%  of  the  freestream  flow  into  the  TWC1  orifice 
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except  that  the  target  plane  (Fig.  2)  was  constrained  to  lie  parallel 
to  the  orifice  plane  at  a distance  1.6  ran  upstream  from  it.  Particle 
flux  tubes  were  centered  at  the  orifice  center. 

To  explore  the  effect  of  flux  tube  radius  on  water  drop 
concentration  factor  results,  calculations  were  done  for  target 
plane  flux  tube  radii  of  0.5,  0.8  and  0.9  inch.  (The  orifice  radius 
is  1 inch.)  Negligible  differences  in  results  were  obtained  over  a 
drop  diameter  range  of  50  to  300  urn.  All  results  reported  below  are 
for  target  plane  flux  tube  radii  of  0.9  inch. 

Concentration  factors,  CF,  for  water  drops  are  given  in 

★ 

Table  7.  Figure  10  shows  a flux  tube  of  10  um  diameter  waterdrops 
to  the  TWCI.  As  seen  in  the  figure,  the  tube  flares  open  as  it  ap- 
proaches the  TWCI;  therefore  we  have  A A,  and  Cp  < 1 (see  eq.  (3)). 
Figure  11  shows  that  otherwise  there  is  little  flux  tube  distortion. 

Concentration  factor  does  not  provide  an  approximation  to 
particle  concentration  distortion  in  the  same  way  that  it  does  for 
the  more  conventional  cases  involving  sampling  instruments  protruding 
from  the  fuselage  aft  of  the  nose.  This  is  because  the  TWCI  orifice 
is  very  nearly  a stagnation  area,  and  we  cannot  assume,  as  usual, 
that  the  particle  velocity  is  essentially  the  same  as  the  air  velocity. 
The  ratio  of  particle  concentration  at  the  TWCI  orifice  to  the  free- 
stream  concentration  is 


C 


M 


= C 


F 


(5) 


where  V . is  the  particle  speed  into  the  orifice  and  V is  freest^eam 
p,t  K 

airspeed.  Table  7 also  contains  the  ratio  V/Vp  t and  C^.  Except 


* Water  drop  concentration  factors  also  were  calculated  for  zero 
intake  into  the  TWCI.  Cp  values  are  not  significantly  different 

from  those  given  in  Table  7. 
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for  the  smallest  water  drops  there  is  little  concentration  distortion 
at  the  orifice.  Obviously,  this  conies  about  because  flux  decrease  is 
compensated  by  concentration  buildup  in  the  staqnation  reoion. 

Though  particle  concentrations  at  the  TWCI  orifice  are 
nearly  the  same  as  freestream  concentrations  this  does  not  mean  that 
hydrometeor  concentrations  measured  by  the  TWCI  are  free  of  distor- 
tion. The  TWCI  measures  the  throughput  rate  of  bulk  particle  mass, 
m . This  is 


= C.A  V . 
t o p,t 


(6) 


where  C^.  is  the  actual  concentration  at  the  orifice,  and  A is  orifice 
area.  From  the  definition  of  concentration  ratio,  the  freestream  con- 
centration, C,  is 


C = 


CMA  V . 
M o p,t 


and  substitution  of  eq.  (5)  into  this  gives 


C 


_m 1 

A V Cr 
o F 


(7) 


Thus  for  the  TWCI,  distortion  in  hydrometeor  concentration  measure- 
ments is  proportional  to  the  reciprocal  of  concentration  factor,  and 
equation  (7)  should  be  used  to  compute  freestream  particle  concentra- 
tions from  the  TWCI  data. 

Concentration  factor  calculations  also  were  done  for  solid 
and  hollow  ice  columns,  hexagonal  plates,  plane  dendrites  and  aggre- 
gates of  unrimed  radiating  assemblages  of  plates,  side  planes,  bullets 
and  columns.  Properties  of  the  ice  hydrometeors  are  given  in  Tables 
4, 5, 8, 9.  (Explanations  of  these  tables  are  in  ref.  2)  Results  are 
given  in  Tables  10-13  and  Fig.  12. 


i 
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TABLE  8 


PROPERTIES  OF  PLANE  DENDRITIC  (Pie)  ICE  CRYSTALS 


Dendri te 
Diameter 
(um) 

Thi ckness 

(um) 

Mass 

U) 

Diameter  Of 
Water  Drop 
Of  Equivalent 
Mass  (um) 

Equivalent 

Densi ty 
(kg/m3) 

Disc 

Diameter 

(um) 

500 

29.16 

0.7356 

112 

334 

341.0 

600 

31.24 

1.0981 

128 

357 

389.4 

800 

34.82 

2.1850 

161 

402 

490.3 

1000 

37.87 

3.6484 

191 

450 

574.1 

1500 

44.13 

8.2800 

251 

556 

720.8 

1800 

47.27 

11.494 

280 

607 

785.4 

2500 

53.50  * 

21.129 

* 343 

~ 712 

924.1 

TABLE  9 

PROPERTIES 

OF  AGGREGATES 

OF  UNRIMED 

RADIATING  ASSEMBLAGES  OF 

PLATES,  SIDE  PLANES,  BULLETS  AND  COLUMNS 


mension 

(pm) 

Mass 

(pg) 

Densi ty 
(kq/m3) 

Diameter  of 
Water  Drop  of 
Equivalent  Mass 
(pm) 

300 

3.756 

265.7 

193 

500 

9.914 

151.5 

266 

800 

24.21 

90.3 

359 

1000 

37.0 

70.7 

413 

2000 

138.1 

33.0 

641 

3000 

298.3 

21.1 

829 
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TABLE  11 


HEXAGONAL  PLATE  ( PI  a ) CONCENTRATION  RATIOS  AT  THE  LEARJET  36  TWCI 


Plate  * 
Diameter 
(um) 

Jr. 

** 

v/vP,t 

CM 

C ” 1 

10 

.401 

3.448 

1.383 

2.492 

30 

.533 

2.433 

1.299 

1.875 

50 

.560 

2.320 

1.297 

1.787 

100 

.601 

2.137 

1.283 

1.665 

300 

.761 

1.558 

1.185 

1.314 

500 

.836 

1.348 

1.128 

1.196 

1000 

.919 

1.151 

1.057 

1.088 

See  Table  5 for  hexagonal  plate  properties. 
See  footnote  in  Table  7. 
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TABLE  12 


PLANE  DENDRITE  (Pie) 

CONCENTRATION  RATIOS 

AT  THE  LEARJET 

36  TWCI 

Dendrite  * 
Diameter 
(-m) 

Jr 

★ ★ 

V / V . 

_P_  ?_t 

C -1 

lf 

500 

.789 

1.472 

1.161 

1.267 

600 

.810 

1.413 

1.145 

1.234 

800 

.842 

1.330 

1. 120 

1.187 

1000 

.866 

1.273 

1. 102 

1.155 

1500 

.900 

1.  190 

1.071 

1.111 

1800 

.912 

1.164 

1.062 

1.096 

2500 

.930 

1.125 

1.046 

1.076 

* 

Toe  Table  8 for 

★ ★ 

See  footnote  in 

plane 

Table 

dendrite  crystal 
7. 

properties. 

TABLE  13 

CRYSTAL  AGGREGATE 

CONCENTRATION  RATIOS  AT  THE  LEARJET 

36  TWCI 

★ 

Aggregate 

Dimension 

(um) 

Jr. 

★ * 

V/Vn  t 

P » l 

CM 

C _1 
_F_ 

300 

.951 

1.058 

1.006 

1.052 

500 

.959 

1.046 

1.002 

1.043 

1000 

.967 

1.033 

.999 

1.034 

3000 

.975 

1.021 

.996 

1.026 

★ 

See  Table  9 for  crystal  aggregate  properties. 

★ ★ 

See  footnote  in  Table  7. 
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Figure  12.  Reciprocal  concentration  factor  vs.  particle  ass  tor  v-r  ous  hydrnmeteors  a 
collected  by  the  Learjet  36  TWCI.  Tables  4, 5, 8, 9 list  nartiule  dimensi'rs  as  • ‘‘^n-.tion 
of  mass. 


mmmm 


CM  .aloes  indicate  teat  ice  parti:  le  concentrations  at 

n 

the  cr;fice  also  are  not  -’uch  different  >o  freestrem  values. 


However , with  exceptions  or  the  larnest  hydro  eteors  of  ail  tyoes 
and  ice  crystal  aqqrenates  of  all  sizes,  Cr  values  ire  substan- 
stially  less  than  one.  This  neans  that  corrections  anr;lied  to  the 
TWCI  data  via  the  reciprocal  concentration  factor  in  eq.  (7) 
usually  are  very  sionificant. 
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GENERALIZATION  OF  THE  CONCENTRATION  FACTOR 
METHOD  FOR  FUSELAGE-MOUNTED  INSTRUMENTS 


INTRODUCTION 

The  concentration  factor  calculation  method  described  above 
and  in  refs.l  and  2 is  designed  to  be  specific  with  respect  to  air- 
plane, sampling  site  and  flight  conditions.  Its  application  requires 
development  of  a detailed  digital  description  of  the  airplane  fuse- 
lage followed  by  a rather  demanding  computational  effort  on  a large 
computer.  Frequently  it  is  desirable  to  make  a preliminary  estimate 
of  the  suitability  of  a sampling  location  without  resorting  to  this 
detailed  analysis.  To  accomodate  such  need,  the  analysis  described 
below  was  done.  It  results  in  a graph  that  can  be  used  to 
roughly  bracket  the  edge  of  the  shadow  zone  for  sampling  points  near 
any  fuselage. 

CONCENTRATION  FACTOR  RESULTS  FOR  POTENTIAL  FLOW  AROUND  AN  ISOLATED 
POINT  SOURCE 

A single  point  source  in  a potential  flow  field  affords 
a simple,  easily  computed  approximation  to  flow  around  an  airplane 
fuselage.  Thus,  this  approximation  was  used  for  the  extensive  con- 
centration factor  calculations  required  for  the  generalization. 

Figs.  13-15  show  water  drop  trajectories  around  two  point  source 
bodies.  The  smallest  drops  (10  ym  diameter)  follow  the  streamflow 
very  closely  and  Fig.  13  essentially  shows  air  streamlines.  Fig.  14 
illustrates  the  increase  with  distance  aft  of  the  depth  of  the 
shadow  zone  for  100  ym  diameter  drops,  and  Fig.  15  shows  the  rela- 
tively unaffected  trajectories  of  large  drops.  The  point  source 
body  definition  and  the  flow  field  equations  are  given  in  the  Appen- 
dix. 
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Trajectories  of  10  diameter  water  drops  in  potential 
airflow  around  a single  point  source  which  simulates  a fuselaqe  of 
radius  0.8  m with  velocity  130  n sec*’  at  20  kft  altitude  Thlse 
trajectories  virtually  coincide  with  airflow  streamlines. 
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a ^Dotentiai  °V°°  diameter  water  drops  through 

a potential  flow  field  around  a single  point  source.  These  cal- 
ulations  simulate  a fuselage  of  radius  0.8  m with  velocity 
130  m sec’1  at  20  kft  altitude. 


Figure  15.  Trajectories  of  1000  pm  diameter  water  drops  in 
potential  airflow  around  a single  point  source.  These  calcula- 
tions simulate  a fuselage  of  radius  0.8  m with  velocity  130  m sec 
at  20  kft  altitude.  The  tilt  of  the  trajectories  is  caused  by 
gravity  settl ing. 


Concentration  factors  were  calculated  such  as  to  bracket 
the  location  of  the  shadow  zone  edge  at  nine  x/D  values  ranging 
from  0.25  to  4.5.  Here  x is  distance  aft  from  the  nose  tip  and  D 
is  the  body  diameter  at  x = ®.  Calculation  geometry  in  the  plane 
normal  to  the  fuselage  axis  is  shown  in  Fig.  16.  The  full  three- 
dimensional  concentration  factor  code  for  water  drops  was  used. 

The  calculations  were  done  for  a standard-atmosphere  alti- 
tude of  20  kft  (T  = 248.6  K,  p = 0.653  kg  m-3),  true  airspeed  of 
130  m sec-1  and  two  fuselage  radii:  0.8  m and  2.0  m.  The  smaller 
fuselage  radius  represents  the  smaller,  more  slender  airplanes  such 
as  the  Cessna  Citation  and  Learjet.  The  larger  radius  represents 
large  airplanes  such  as  the  Lockheed  C130.  The  only  significant 
difference  between  results  for  the  two  fuselages  is  that  the  concentra- 
tion factor  vs.  drop  diameter  curves  peak  at  about  50  urn  for  the 
fuselage  of  smaller  radius  while  they  peak  at  about  80  pm  for  the 
larger  fuselage. 

It  is  not  possible  to  precisely  determine  the  shadow  zone 
edge  via  this  method.  This  is  because  as  the  shadow  zone  is  ap- 
proached concentration  factors  increase  sharply  and  numerical  prob- 
lems ensue.  Thus,  it  is  only  possible  to  rather  broadly  bracket  the 
shadow  zone  edge.  The  stippled  region  of  Fig.  17  represents  the 
bracketed  region:  r-r  values  (r-r  corresponds  to  y-yc  in  Fig.  16) 
above  the  stippled  area  are  clear  of  the  shadow  zone,  values  below 
the  stippled  area  are  within  the  shadow  zone,  and  values  within  the 
stippled  area  may  or  may  not  be  within  the  shadow  zone  but  if  not 
will  have  extremely  high  concentration  distortion.  The  graph  brack- 
ets the  maxima  of  the  concentration  factor  vs.  water  drop  diameter 
curves.  Distances  are  normalized  by  dividing  by  maximum  body 
diameter,  D. 
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Figure  16.  Cross  section  of  simulated  fuselage  showing  concen- 
tration factor  calculation  geometry  used  for  the  generalization 

calculations.  SP  is  the  sampling  point,  y-y  is  the  distance 

s 

of  the  sampling  point  from  the  fuselage  surface  measured  along 
the  horizontal  (y)  axis.  The  z axis  is  vertical. 


L 


Figure  17.  Summary  of  concentration  factor  results  that  bracket  the  shadow 
zone  edge,  x is  distance  aft  from  the  nose  tip. 

fuselage  surface.  (For  the  flow  around  the  point  source  body. 


r-rs  is  distance  from  the 


3 y-y< 


as  in  Fig.  16.)  D is  maximum  fuselage  diameter.  The  stippled  area  brackets 
the  shadow  zone  edge  for  flow  around  the  point  source  body.  Case  study 
points  are: 


0 

C130A 

Formvar  Replicator 

* 

Cessna  Citation 

0 

C130E 

Formvar  Replicator 

A 

C130B 

tyman-n  Total  Water 

• 

Learjet  36 

V 

Probe 

■ 

Lear jet  36 

C130B 

Foil  Impactor 

A 

Learjet  36 

PMS  Precipitation 
Particle  Spectrometer 


PMS  Cloud  Particle 
Spectrometer 

Snow  Stick 
Formvar  Replicator 

PMS  Axial  Scattering 
Probe 


See  text  for  additional  description. 
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COMPARISON  WITH  CASE  STUDY  RESULTS 


Table  14  summarizes  pertinent  results  of  the  case  studies 
performed  to  date.  Only  for  the  Cessna  Citation  are  the  results  for 
the  instrument  sampling  points;  for  the  other  cases  the  results  are 
for  points  closer  to  the  fuselages  along  the  instruments'  axes.  Where 
available,  results  that  bracket  the  shadow  zone  edge  are  qiven. 

(For  the  C130A  formvar  replicator  and  Cessna  Citation  PMS  probes, 
shadow  zone  edges  were  not  bracketed  and  results  for  single,  clos- 
est points  are  given.)  These  data  also  are  plotted  in  Fig.  17. 

Plotted  points  that  bracket  the  snadow  zone  are  connected  by  ver- 
tical lines  in  the  figure.  The  number  above  the  upper  point  of  a 
pair,  which  is  clear  of  the  shadow  zone,  is  the  largest  concentra- 
tion factor  calculated. 

Figure  18  shows  normalized  fuselage  radii  values  as  a 
function  of  normalized  distances  from  the  nose  tip  for  the  various 
airplanes  and  for  the  point  source  body.  For  the  airplanes,  the  y$ 
values  are  for  plan-view  projections  of  the  fuselages.  The  figure 
shows  that  the  point  source  body  is  the  bluffest  obstruction  to 
flow,  but  affords  a reasonable  approximation  to  the  fuselage  shapes 
with  the  exception  of  the  Learjet. 

That  the  point  source  body  has  the  bluffest  nose  shape  may 
at  least  partly  explain  why  the  case-study  data  all  fall  below  that 
for  the  point  source  body  for  nose-mounted  instruments  (x/D  < 2). 

The  generalized  and  case-study  data  for  instruments  mounted  on  the 
sides  of  the  cylindrical  sections  of  fuselages  (x/D  > 2)  agree  very 
wel  1 . 

Figure  17  should  be  considered  as  a preliminary  result; 
more  exhaustive  calculations  are  necessary  before  the  generaliza- 
tion can  be  considered  to  be  complete.  Nevertheless,  the  figure 
may  be  used  to  provide  rough  estimates  of  "safe"  water  drop  sam- 
pling distances  from  fuselage  surfaces. 
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TABLE  14 

SUMMARY  OF  CASE-STUDY  CONCENTRATION  FACTOR  RESULTS 


Airplane 
(Fuselage 
Diameter  (m)) 

Instrument 

x/D 

r-r 

s 

__D 

★ 

r 

'F 

C130A 

(4.32) 

Formvar 
Repl icator 
(Ref.  1) 

0.682 

0.061 

1.7 

C130E 

(4.32) 

Formvar 
Repl icator 
(Ref.  1) 

1.627 

0.104 

0.081 

2.5 

S 

C130B 

Lyman— i 

1.277 

0.075 

1.8 

(4.32) 

Total  1-Jater 
Probe  (Ref.  11) 

0.057 

S 

C130B 

Foil  Impactor 

1.562 

0.085 

2.0 

(4.32) 

(Ref.  11) 

0.063 

S 

Cessna 

Citation 

(1.61) 

PMS  Precipitation 
Particle  Spectrometer 
(Ref.  1) 

2.632 

0. 142 

1.6 

Cessna 

Citation 

(1.61) 

PMS  Cloud 

Particle  Spectrometer 
(Ref.  1) 

2.632 

0.142 

s 

Learjet  36 
(1.62) 

Snow  Stick 

3.133 

0. 164 
0.124 

1.4 

$ 

Learjet  36 

Formvar 

2.720 

0.164 

1.3 

(1.62) 

Repl icator 

0.124 

s 

Learjet  36 

PMS  Axial  Scattering 

2.384 

0.140 

1.3 

(1.62) 

Probe 

0. 100 

s 

★ 

Concentration  factors  are  peak  values  for  water  drops. 

+0nly  for  the  Cessna  Citation  PMS  probes  are  the  results  given  for  the 
instrument  sampling  points.  In  all  other  cases  the  concentration  factor 
results  are  for  points  closer  to  the  fuselages  along  the  instruments'  axes. 

★ 

S denotes  that  this  sampling  point  lies  inside  the  shadow  zone. 


- 49  - 


Point  Source 

C130A 

C130E 

Cessna  Citation 
Learjet  36 


Figure  13.  Fuselage  radius,  y vs.  distance  aft  from  the  nose  tip 

x,  for  various  airplane  fuselaqes  and  the  point  source  body.  D is 
maximum  fuselage  diameter.  See  text  for  additional  description. 
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Unfortunately  there  are  not  enouqh  ice  hydrometeor  results 
available  to  warrant  construction  of  similar  graphs  for  ice.  How- 
ever, the  reader  should  be  aware  that  shadow  zones  are  deeper,  extend 
to  larger  particle  masses  and  span  broader  ranqes  of  sizes  for  ice  as 
compared  with  water  drops  (see  refs.  1 and  2). 
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SUMMARY  AND  CONCLUSIONS 


Concentration  factors  for  water  drops  were  computed  at  the 
sampling  points  of  three  instruments  mounted  through  the  cabin  win- 
dows of  the  Learjet  36  airplane:  a PMS  axial  scattering  probe,  a snow 
stick  and  a formvar  replicator.  The  largest  concentration  factor  is 
1.2  for  80  urn  diameter  drops  at  the  snow  stick,  which  instrument  is 
the  furthest  aft  of  the  three.  To  determine  if  the  snow  stick  is 
shadowed  for  ice  hydrometeors,  concentration  factors  for  ice  columns 
and  plates  also  were  computed  for  that  instrument;  no  shadowing  was 
indicated.  The  distance  of  the  shadow  zone  edge  from  the  fuselaoe 
surface  was  determined  approximately  for  each  instrument. 

Collection  efficiency  for  hydrometeors  of  the  Total  Water 
Content  Indicator  (TWCI)  mounted  in  the  Learjet  36  nose  was  evaluated. 
To  accomplish  this,  the  concentration  factor  code  was  modified  to 
allow  intake  of  a specified  fraction  of  the  freestrean  flow  by  the 
TWCI  orifice.  Detailed  examination  of  the  air  flowfield  into  and 
around  the  orifice  showed  that  the  computed  flow  is  highly  accurate. 
Hydrometeor  flux  tubes  were  calculated  into  the  TWCI  orifice  for  water 
drops  and  various  ice  crystal  forms.  Concentration  factors  substan- 
tially less  than  unity  are  found  for  small  hydrometeors  of  all  types 
except  for  ice  crystal  aggregates.  Since  freestream  hydrometeor 
concentrations  determined  by  the  TWCI  are  proportional  to  reciprocal 
concentration  factor  (eq.  (7)),  larqe  correction  factors  are  reauired 
for  most  small  hydrometeors  (Fig.  12). 

A preliminary  version  of  a generalized  graph  by  which  "safe" 
water  drop  samplino  distances  for  fuselage-mounted  instruments  can  be 
estimated  was  prepared  (Fig.  17).  Potential  airflow  around  an  ideal 
point  source  was  used  as  a convenient  model  of  an  airplane  fuselaoe 
for  the  extensive  concentration  factor  calculations  required.  Re- 
sults for  the  point  source  body  are  compared  with  available  case 
study  results.  Agreement  is  good  for  instruments  mounted  on  the 
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cylindrical  portions  of  fuselages.  For  instruments  mounted  on  the 
forward  (tapering)  parts  of  fuselages,  case  study  results  indicate 
more  shallow  shadow  zones  than  do  the  point  source  body  results. 
This  latter  may  be  caused,  at  least  in  part,  by  the  fact  that  the 
nose  of  the  point  source  body  affords  the  bluffest  obstacle  to 
flow  of  all  shapes  considered. 


APPENDIX 


VELOCITY  FIELD  EQUATIONS  FOR  POTENTIAL  FLOW 
AROUND  A SINGLE  POINT  SOURCE 


The  following  is  developed  from  chapter  16  of  the  book  by 
Mil  ne-Thompson^^ . The  stream  function,  iji,  for  unit  freestream 
flow  about  a source  of  strength  1/4  at  the  origin  is 


4q 


(A.  1 ) 


where 


q2  = x2  + y2  + z2 


r2  = y-  + z2 


and  freestream  flow  is  parallel  with  the  x axis  and  toward  the  posi- 
tive x direction. 

The  velocity  equations  are 


(A. 2) 


1 


4q 


The  body  surface  (defined  by  = - 


s given  by 


(A. 3) 


x 

^s 


) 


12.  L.M.  Milne-Thompson,  Theoretical  Hydrodynamics,  fifth  edition 
(The  MacMillan  Co.,  1968T- 
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which  can  be  solved  to  give  1%  as  sinole  valued  function  of  x: 

V y 2 (x  vx'  + 2 + 1 - X ) ; X - - I . (A. 4) 

Equation  (A. 4)  is  used  to  calculate  the  surface  contour  of  the  body. 
The  nose  or  stagnation  point  is  at  x = - ],  and  the  body  has  unit 
radius  (r$  = 1)  at  x = ». 
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